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Abstract 

The pseudorapidity dependence of elliptic flow in Au+Au collisions at 130 
A GeV is studied within a full three-dimensional hydrodynamic model in the 
light-cone coordinate. First, we prepare two initial conditions in the hydro- 
dynamic model for analyzing elliptic flow. Both initial conditions lead to 
reasonable agreement with single particle spectra in central and semi-central 
collisions. Next, by using these hydrodynamic simulations, we compare el- 
liptic flow as a function of pseudorapidity with experimental data recently 
measured by the PHOBOS Collaboration. Our results are in agreement with 
experimental data only near midrapidity. This suggests that thermalization in 
the early stage of collisions is not achieved in forward and backward rapidity 
regions. 
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Nucleus-nucleus collisions at the Relativistic Heavy Ion Collider (RHIC) give us an op- 
portunity to study a new state of deconfined nuclear matter, the quark-gluon plasma (QGP) 
|JI|]. The main goals in the physics of relativistic heavy-ion collisions are not only the dis- 
covery of the QGP, but also the investigation of thermodynamical aspects of its new phase, 
i.e., the equation of state (EOS), the order of phase transition between the QGP phase and 
the hadron phase, or the critical temperature. It is very important to check whether the 
thermalization of the nuclear matter produced at the very early stage of collisions is really 
achieved, before discussing the thermodynamics of the QGP. Elliptic flow in non-central 
collisions |^ is suited for this purpose. If produced particles frequently rescatter with each 
other, we naively expect that thermalization of the system is achieved and that the large 
pressure is built in the reaction zone. The pressure produces the momentum anisotropy 
of observed particles from the spatial deformation in the transverse plane. The coefficients 
of second harmonics in the azimuthal distribution, not only its magnitude V2 but also its 
transverse momentum dependence V2{pt) and its centrality dependence V2{h), seem to be 
good indicators for thermalization of nuclear matter Hydrodynamic simulations [0, 

in which full local thermalization is assumed, give us excellent agreement with the first re- 
sult of V2{pt) from Au-|-Au 130 A GeV collisions at the RHIC up to ~ 1.5 GeV near 
midrapidity. Their results are based on a (2+l)-dimensional hydrodynamic model with the 
Bjorken's scaling solution. One cannot, however, discuss the rapidity dependence of observ- 
ables by exploiting the Bjorken's model |Tl|. Since this model is assumed to give a good 



description of space-time evolution of nuclear matter near the midrapidity region, the agree- 
ment between the model calculation and experimental data means that early thermalization 
at the RHIC is achieved at least only at midrapidity. In this paper, we investigate how 
far from midrapidity thermalization in the early stage is achieved through comparison the 
pseudorapidity dependence of elliptic fiow from a genuine three-dimensional hydrodynamic 
model with experimental data recently measured by the PHOBOS Collaboration ||12|| . 

There are several works based on full three-dimensional hydrodynamic simulations for 
one fiuid [1T3|-[T7[]. These simulations, except Ref. 



[|17| , are performed in the Cartesian coordi- 
nate. In view of numerical analysis, hydrodynamic simulations in the Cartesian coordinate 
at the RHIC energy are tremendously hard owing to very long life time of the fiuid (t/ ~ 100 
fm). So we extend our previous algorithm in the Cartesian coordinate [|I6| to the 'light-cone' 
(or Bjorken) coordinate for saving of the computation time. In the new coordinate, hydro- 



dynamic equations [d^T^^^ = 


D and c^^n^ = 0) are written as 
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Here E, P, and ub are energy density, pressure, and baryon density; r = y/2x+x^ = 
V^P"^^^, rjs = (1/2) log(x+/x_) = {l/2)\og[(t + z)/ (t — z)], x, and y are, respectively, 
the proper time, the space-time rapidity, the transverse coordinate parallel to the impact 
parameter vector, and the coordinate perpendicular to the reaction plane. Note that V = 
{d^, dy, driJr). Fluid velocities in the new coordinate are represented in terms of v^, Vy, and 
Vz', Vx = cosh If/ cosh(lf — r/s), Vy = f^coshYf/ cosh(yf — rjs), and v^j^ = tanh(Yf — ?7s), 
where Yf = (1/2) log[{l+Vz)/{l — Vz)] is the rapidity of a fluid element. The Lorentz gamma 



factor is 7 = — v^. We use the same model EOS as the one represented in Ref. |]T5| 



The EOS has a first order phase transition between the QGP phase and the hadron phase 
at Tc(?T,B = 0) = 160 MeV. The QGP phase is assumed to be free gas composed of quarks 
with Nf = 3 and gluons. For the hadron phase we adopt a resonance gas model, which 
includes all baryons and mesons up to the mass of 2 GeV []18[, together with an exclude 



volume correction |jT9|. For further details on the EOS, see Refs. [|n 

We set up initial conditions at tq = 0.6 fm. The initial energy density is motivated by a 
tilted disc recently discussed in the context of 'anti-flow' or 'third flow component' ||20| , |2T| . 
We assume that the longitudinal proflle of initial energy density E{7]s) at a transverse coor- 
dinate (x, y) is composed of two regions; the initial energy density is flat near ~ and 
smoothly connects to vacuum as a half part of a Gaussian function in the forward and back- 
ward space-time rapidity regions. The length of a flat region Ar^flat and the width of Gaussian 
function Ar^causs are adjustable parameters to be determined by the experimental data of 
the (pseudo)rapidity distribution. In symmetric collisions with zero impact parameter, we 
expect E{ris) = E{—ris). On the other hand, we shift the energy density by Arjs which is 
identifled with the center-of-rapidity for each transverse coordinate |2^ in non-central colli- 
sions. The difference of thickness between two colliding nuclei for each transverse coordinate 
results in the initial distribution like a tilted disc. We compare two initial conditions in this 
paper. In order to reproduce a dip structure of the pseudorapidity distribution observed by 
the PHOBOS Collaboration we choose A%at = 6.0, Ar/causs = 0.4, and a maximum 
energy density -Emax = 40 GeV/fm'^ for the initial condition A (IC A). The other initial con- 
dition (IC B) has a smaller flat region Ar^flat = 2.8, a larger Gaussian width A?7Gauss = 1-6, 
and a slightly larger maximum energy density E^^^ = 43 GeV/fm^ than the IC A. The latter 
also gives a reasonable result of the pseudorapidity distribution. Both initial conditions for 
energy density at t^s = are the same as the model 'eBC in Ref. E^, i.e., the transverse 



proflle of energy density is in proportion to the number of binary collisions in the transverse 
plane. This model gives the similar shape of the charged particle yield per participating 
nucleon pair ((iA'ch/'i'7)/(0.5A^part) as a function of A^part O]- Figure 1 shows two sets of 



initial energy density in the reaction plane. We adopt the following initial baryon density 
and flow velocities as common initial conditions. For the initial baryon density, we extend 
the parameterization by Sollfrank et al. ||22|, originally proposed for central collisions, to the 
case for non-central collisions. We found that the shape of baryon density does not largely 
affect the final results obtained in this paper since the produced matter at the RHIC energy 
is dominated by mesons. Initial longitudinal flow is the Bjorken's solution Vr,X'To) = (or 
Yf(ro) = rjs) and transverse flow velocities vanish VxiTo) = Vy^To) = 0. 

Figure 2 shows the pseudorapidity distribution of charged particles for central and semi- 
central collisions. We accumulate the contribution from particles directly emitted from 
freeze-out hypersurface and feeding from resonance decays. The preliminary experimental 
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data are given by the PHOBOS Collaboration [^. The PHOBOS Collaboration estimates 
the average numbers of participants < A^part >~ 340 for the 0-6% central collisions and 



93 for the 35-45% semi-central collisions EM. We choose the impact parameters 2.4 fm 



for central and 8.9 fm for semi-central collisions. The resultant numbers of participants based 
on the wounded nucleon model with the standard Woods-Saxon profile are 342 {b = 2.4 fm) 
and 94 {b = 8.9 fm), respectively. We first tune initial parameters to reproduce central 
events. By changing the impact parameter and unchanging the other parameters, we also 
obtain the results for another centrality. Although the IC A is the best fit to experimental 
data, the IC B is also acceptable within error bars. It should be noted that a dip structure 
at 1] = results from the Jacobian of the transformation from rapidity Y to pseudorapidity 

V 111 



In order to describe the dominant radial flow at freeze-out, we choose an appropriate 
value for freeze-out energy density so as to reproduce the slope of transverse momentum 
distributions. In Fig. 3 we compare pt spectra of negative pions and anti-protons for central 
and semi-central collisions with experimental data measured by the PHENIX Collaboration 
p6| . By choosing the freeze-out energy density Ef = 120 MeV/fm^ which corresponds to the 
mean freeze-out temperature < Tf >|.^^=o ~ -'-^^ MeV, our results are in very good agreement 
with the slope of these spectra. The resultant mean radial flow < "^1. >|^ are 0.49 {b = 2.4 
fm) and 0.41 (6 = 8.9 fm) for the IC A and 0.50 {b = 2.4 fm) and 0.42 \b = 8.9 fm) for the 
IC B. 

Finally we perform hydrodynamic simulations with various impact parameters and obtain 
the pseudorapidity dependence of elliptic flow V2{ri) for charged particles in minimum bias 
events 
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Our results and experimental data obtained by the PHOBOS Collaboration are shown in 
Fig. ^. For the IC A, two bumps appear at \ri\ ~ 3. This structure is caused by the tilted 
disc shape, i.e., produced matter is highly deformed in forward and backward space-time 
rapidity regions (see also Fig. 1 (a)). The initial energy density for the IC A has a 'crescent' 



shape rather than an almond shape at \r]s\ ~ 3.5 in the transverse plane |]2l[. Since the initial 
deformation of energy density for the IC B is almost independent of r]s owing to a relatively 
large A?7Gauss (see also Fig. 1 (b)), the resultant ^2(^7) has no bumps. Thus we find the 
pseudorapidity dependence of elliptic flow is highly sensitive to the initial longitudinal profile 
of energy density. In both cases, our results are in reasonable agreement with experimental 
data only near the midrapidity region. These results indicate that thermalization is fully 
achieved near the midrapidity region and partially in forward and backward regions. We 
note that our V2{pt) and f2(^ch/'^max) in I'?! < 1-3 (not shown) are consistent with the results 
shown in Ref. ||^; our results are consistent with the STAR data |TD| below ~ 1 GeV and 
above nch/n^ax ~ 0.5. 

Three discussions are in order here: 

(i) It is best to compare elliptic flow from hydrodynamics with experimental data for each 
centrality in order to see how the thermalized region increases with centrality. In addition to 
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elliptic flow, we need to analyze the transverse momentum (mass) distribution in the forward 
rapidity region since radial flow is also appropriate for checking thermalization of the system. 
There are no such data at present. Therefore we hope the forthcoming experimental data 
of V2{'rj) for each centrality or the transverse momentum distribution in the forward rapidity 
region will reveal the dynamics of nuclear matter in the whole phase space. 

(ii) It is assumed in conventional hydrodynamic models that the sharp transition from 
the hydrodynamic piture to the non-interacting particle picture happens at the constant 
temperature. On the other hand, the sophisticated freeze-out procedure yields small values 
of V2 even if thermalization sets in the early stage of collisions The model evolves the 
QGP and mixed phase as a relativistic fluid, while it switches to a hadronic cascade model 
(RQMD) at a temperature Tswitch = 160 MeV. The main advantage of this hybrid model is 
that the freeze-out process is automatically included without any assumption such as freeze- 
out temperature. This hybrid model can give small V2 in forward or backward rapidity 
regions where the multiplicity or the number of hadronic rescattering is small. 

(iii) There may exist other initial conditions for hydrodynamic simulations which lead 
to be in agreement with single particle spectra and elliptic flow. Along the lines of thought, 
more sophisticated initialization based on the nuclear thickness function will be studied 
elsewhere P7[] . 

In summary, we have studied early themalization of produced particles at the RHIC by 
using a genuine three-dimensional hydrodynamic model. Our hydrodynamic simulations are 
performed in the light-cone coordinate. Such simulations become important in analyzing 
ultra-relativistic heavy-ion collisions at the RHIC energies. We considered two sets of initial 
parameters in the hydrodynamic model. We chose initial parameters so as to reproduce 
the pseudorapidity and transverse momentum spectra for both central and semi-central 
collisions. By using results from these numerical simulations, we analyzed the pseudorapidity 
dependence of elliptic flow. We found that a large value of V2 comes from the tilted disc 
shape which has a large spatial deformation in the transverse plane in the forward rapidity 
region. Near rj = our results for minimum bias events are consistent with experimental 
data observed by the PHOBOS Collaboration. On the other hand, we failed to reproduce 
experimental data in forward and backward rapidity regions (|?7|^1). This indicates the early 
pressure is not built in forward and backward regions. Both theoretical and experimental 
detailed analyses as discussed above are indispensable for understanding the dynamics of 
nuclear matter not only at midrapidity but also at forward (backward) rapidity at the RHIC 
energy. 

The author is much indebted to S. Esumi, T. Hatsuda, P. Huovinen, Y. Miake, K. Morita, 
S. Muroya, and C. Nonaka for fruitful discussions. He also thanks B. Back, S. Manly, and 
I. Park for teaching experimental aspects to him and C. Nonaka for providing him with a 
numerical table of EOS. He also acknowledges Inoue Foundation for Science for the flnancial 
support. 
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FIGURES 




FIG. 1. Two sets of initial energy density in the reaction plane (y = and % > 0). (a) Initial 
condition A. (b) Initial condition B. See text for details. 
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5 300 




FIG. 2. Pseudorapidity distribution of charged particles in Au+Au 130 A GeV central and 
semi-central collisions. Solid lines and dashed lines correspond to initial conditions A and B, 
respectively. Experimental data are observed by the PHOBOS Collaboration. 
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FIG. 3. Scaled transverse momentum distribution of negative pions and anti-protons in Au+Au 
130 A GeV central and semi-central collisions. Solid lines and dashed lines correspond to initial 
conditions A and i?, respectively. Experimental data are observed by the PHENIX Collaboration. 
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preliminary PHOBOS data 
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FIG. 4. Pseudorapidity dependence of elliptic flow for charged particles in Au+Au 130 A GeV 
collisions. The value of elliptic flow is averaged over all centrality. Solid line and dashed line 
correspond to initial conditions A and B, respectively. Experimental data is observed by the 
PHOBOS Collaboration. 
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